Accepted by ApJ: 25 October 2007 

XMM-Newton Detection of a Compton-thick AGN in the 1-Jy 

ULIRG/LINER F04103-2838 

■|^-, Stacy H. Teng| S. Veilleux, A. S. Wilson 

Q ' Department of Astronomy, University of Maryland, College Park, MD 20742 

-)— • i A. J. Youne 

s-^ \ Center for Space Research, Massachusetts Institute of Technology, Cambridge, MA 02139 

m ; D. B. Sanders 

r-| ' Institute for Astronomy, University of Hawaii, 2680 Woodlawn Drive, Honolulu, HI 96822 

in' 

to N. M. Nagar 

\ Astronomy Croup, Departamento de Fisica, Universidad de Concepcion, Casilla 160-C, 

^ \ Concepcion, Chile 

§ ; ABSTRACT 

p; 

^-H ' We report on the detection of Fe Ka emission in F04103— 2838, an ultralu- 

f~^ ■ minous infrared galaxy (ULIRG; logfLia/Lo] > 12) that is optically classified as 

O ■ a LINER. Previous Chandra observations suggested the presence of both a star- 

^^ ' burst and an AGN in this source. A deeper (~20 ksec) XMM-Newton spectrum 

^ ■ reveals an Fe Ka line at rest frame energy ~6.4 keV, consistent with cold neutral 

c^ ■ iron. The best-fit spectral model indicates the Fe Ka line has an equivalent width 

of ~1.6 keV. The hard X-ray emission is dominated by a Compton-thick AGN 

with intrinsic 0.2-10 keV luminosity ~10^^ ergs s~^, while the soft X-ray emission 

is from ~0.1 keV gas attributed to the starburst. The X-ray spectrum of this 

source bears a striking resemblance to that of the archetypal luminous infrared 

galaxy NGC 6240 despite differences in merger state and infrared properties. 
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Introduction 



The primary energy source (AGN versus starburst activity) of ultraluminous infrared 
galaxies (ULIRGs; logfLiR/L© > 12]) is still a matter of debate. Optical and infrared 
emission-line spectra suggest that the energy output of most local ULIRGs is dominated 
by starbursts, but the "warm" infrared colors and quasar-like spectra of the more luminous 
ULIRGs indicates that black-hole driven activity plays an increasingly important role in 
these objects (e.g., Veilleux et al. 1995, 1997, 1999a, 1999b; Genzel et al. 1998; Surace & 
Sanders 1999; Tran et al. 2001). The dusty, gas-rich nature of ULIRGs implies, however, 
that observations in energy bands other than radio and X-ray may not always probe the 
true nuclear energy source of these objects. Since the luminosity of ULIRGs in the radio is 
insignificant. X-ray observations remain arguably the best option to solve this energy source 
mystery. 

Unresolved hard X-ray emission is in principle a telltale sign of a dominant AGN. How- 
ever, if a large column density of gas (> 10^^ cm~^) is located in front of the nucleus, then 
directly viewed X-rays from the AGN will be strongly attenuated. For such cases, Fe Ka 
lines with large equivalent widths (> 1 keV) are expected due to scattering off circumnuclear 
material (e.g., Ghisellini et al. 1994; Krolik 1994). Thus, the discovery of such Fe Ka lines 
may be the best evidence for energetically dominant AGNs in highly obscured ULIRGs. 



In re cent years, t 



iree X-ray surveys have added considerably to our knowledge of 



ULIRGs. iPtak et al.l (120031 ) performed a volume-limited {z < 0.045) survey of ULIRGs 



with Chandra. On the basis of their dust temperatur es (25-to-60 /xm f lux ratio) and X- 



ray luminosities, three of the eight ULIRGs sampled by iPtak et al.l (12003) were classified as 



AGN- dominated (Mrk 231, Mrk 273, and F05189-2524). In the same year. iFranceschini et al. 



( 120031 ) published the results of a similar survey with XMM-Newton that focused on the bright- 



est local U LIRGs (only one ULIRG in their sample had z > 0.082). Of the ten ULIRGs 



sampled bv IFranceschini et al.l (l2003f l. three were AGN-dominated (Mrk 231, F19254-7245, 
and F20551-4250) and two had X-ray signatures of both a starburst and an AGN (F20100- 
4156 and F23128-5919). All of the AGN-dominated ULIRGs showed strong Fe K emission 
lines (Maloney & Reynolds 2000; Braito et al. 2003; Ptak et al. 2003; Franceschini et al. 
2003; Braito et al. 2004). These two pioneering surveys proved the viability of using the 
X-ray emission as a diagnostic for AGN activity in ULIRGs. However, these two surveys 
only studied a small set of the nearest and brightest ULIRGs, and therefore were not able 
to draw general conclusions on the issue of the energy source among ULIRGs as a class. 

In an attempt to expand this type of study to a more characteristic sample of ULIRGs, 
our group (Teng et al. 2005, hereafter Paper I) conducted a snapshot (10 ksec/target) survey 
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of 14 ULIRGs from the 1-Jy samplcl. These sources were carefully selected to sample the 
full range of infrared luminosities and infrared colors that characterize the entire class of 
local ULIRGs. All 14 galaxies were detected by Chandra, though most (11/14) had less than 
40 counts. The analysis showed that the two brightest galaxies in the sample have optical 
and X-ray spectral characteristics of Seyfert 1 nuclei. Most others have X-ray photon indices 
(estimated using hardness ratios) and hard X-ray to far-infrared flux ratios which are similar 
to those of starbursts. 

One exception, F04103— 2838, had a hardness ratio (deduced from only 30 counts) that 
suggested the presence of a starburst coexisting with an AGN. The low signal-to-noise data 
could not distinguish between a Compton-thick AGN or an intrinsically faint nuclear source. 
This object is optically classified as a LINER (Veilleux et al. 1999a). F04103-2838 has 
one of the largest 25-to-60 /im flux ratios of all 1 Jy ULIRGs (/25//60 = 0.30). In fact, this 
is the warmest of all IRAS 1 Jy ULIRGs with optical LINER or H II classification. This 
source is even warmer than some of the Seyfert galaxies in the 1 Jy sample (see Fig. 1 of 
Paper I). Recent Spitzer detection of [Ne V] and [O IV] lines from this source has confirmed 
the existence of an AGN in this system (Veilleux et al. 2007, in prep.), making it a rare 
example of a ULIRG optically classified as a LINER that is not classified as a starburst 
on the basis of mid-infrared spectroscopy (Lutz, Veilleux, & Genzel 1999; Taniguchi et al. 
1999). In this paper, we present an XMM-Newton observation of F04103— 2838 which delves 
deeper into the nature of this AGN. 

F04103— 2838 is an interacting galaxy system in the l ate stages of a merger as indicated 



by the presence of a single nucleus with distinct tidal tails ( IVeilleux et al.ll2002l : iDasyra et al. 



2006al ). This object has an infrared (8-1000 //m) luminosity of lO^^'^^L© and a cosmolog- 
ical redshift of 0.118. Assuming Hq = 75 km s^^ Mpc^^ and qo = (used throughout 
this paper), the luminosity distance of this object is 497 Mpc. At this distance, 1" corre- 
sponds to ~2.4 kpc. In §2 of this paper, we describe our new XMM-Newton observation of 
F04103— 2838 and the methods we used to reduce these data. In §3, we present the analysis 
of these data, emphasizing the results on the X-ray morphology, the lack of flux variability, 
and the spectral decomposition of the X-ray emission. The implications of these results are 
discussed in 54. The main conclusions are summarized in 55. 



^The 1 Jy sample of ULIRGs is comprised of IRAS galaxies with fluxes at 60 ^m exceeding 1 Jy, 
LiR > 10^^ L0, galactic latitude |b| > 30°, f(60 fim) > f(12 /im) (to avoid stars), IRAS color log[/6o//ioo] 
> -0.3 (to favor luminous infrared systems), and redshift 0.018 < z < 0.268 (Kim & Sanders, 1998). 
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2. Observation and Data Reduction 

F04103— 2838 was observed with XMM-Newton during orbit #1132 on 13 February, 
2006 (ObsID: 0301330401; PI: Wilson) with the EPIC instrument. The EPIC cameras were 
operating in full-frame mode. Each of the detectors used the medium filter. 

The data were processed using the standard procedures of the XMM-Newton Science 
Analysis System (SAS) version 6.5.0 released on 17 August, 2005. The processing procedures 
outlined in §4.11 of the XMM-Newton SAS User's Guide were followed. The event lists 
were calibrated with the latest available calibration files as of June, 2006. Times of high 
background flares were flagged. The total good time interval on source for each camera was 
17.5 ksec for PN, 21.8 ksec for MOSl, and 21.5 ksec for M0S2. 

Source and background counts were extracted from circular regions with radii of 24". 
Because the source is near a gap in the CCD and a nearby X-ray luminous source, the 
background was extracted from a circular region with the same area as the source extraction 
region in a neighboring piece of the sky in which no obvious X-ray sources reside. The 
total 0.2-10 keV counts extracted from the source region are 224 for PN, 52 for MOSl, and 
48 from M0S2. The expected background counts in the source region are 45 for PN and 
50 for MOSl/2 based on the expected background count rates quoted in the XMM User's 
Handbook. Since the extracted source counts are approximately the same as the expected 
background counts for the MOS detectors, we will exclude the MOS spectra in our spectral 
analysis of the source. 



3. Analysis 

In H'd.l\ we describe the distribution of the X-ray emission from F04103— 2838. In §3.21 
we point out the lack of variability of this object. A detailed analysis of the X-ray spectrum 
and iron complex is presented in §3.3[ 



3.1. Morphology 

To improve the signal-to-noise ratio of the images, the PN and MOS 1/2 events were 
combined using the SAS task emosaic and then smoothed with a 5" Gaussian using asmooth 
to match the spatial resolution of XMM. The resultant image is displayed in the left panel of 
Figure [TJ A comparison of the 0.2-2 keV (unsmoothed) radial profile with the XMM-Newton 
point spread function (PSF) at 1 keV indicates that the source is unresolved (see Fig. 1, 
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right panel). Only the EPIC PN data were used for the radial profile calculations because of 
the small number of counts detected by the MOSl/2 cameras. The PSF of the PN camera is 
well described by a King profilcl and was normalized so that the total number of counts per 
square pixel under the curve match the total number of detected counts per square pixel. 
The Chandra data from Paper I verifies that the source is unresolved. 



3.2. X-Ray Variability 

The time interval covered by our observation was divided into four equal bins of 5234 
seconds to search for significant X-ray variability, another potential indicator of dominant 
AGN activity. The 0.2-10 keV and 2-10 keV EPIC PN count rates were calculated for both 
source and background. Figure [2] shows the 0.2-10 and 2-10 keV light curves of the source 
and background. To within the errors, the source is not significantly variable on the 5-6 hour 
time scale of our observations. 



3.3. X-Ray Spectra 

The extracted source and background spectra from each detector were binned using the 
FTOOL grppha to at least 3, 5, and 15 counts bin~^. The binned and unbinned spectra 
were then analyzed using XSPEC version 11.3.2t. The quoted errors on the derived best- 
fitting model parameters correspond to a 90% confidence level (A^^/Ac-stat = 2.706). The 
X^ goodness-of-fit test was used to judge the fits to the spectrum binned to at least 15 
counts bin~^. The Cash statistics (c-stat) option in XSPEC was used for spectra binned to 
at least 3 and 5 counts bin~^ and the unbinned data. The spectral model was applied to 
the EPIC PN data only (see ||2]). Al l models were corrected f or Galactic absorption using 



NH,Gaiactic = 2.45 X 10^° atoms cm ^ ( iDickey fc Lockmadll990l ) 



3.3.1. Effects of Binning 

By definition, spectra binned to at least 15 counts bin~^ have the highest signal-to-noise 
ratios while the spectra binned to at least 3 counts bin~^ show the most spectral details. 
The first task is to determine whether the mode of binning affects the spectral parameters 



^PSF == A[l + (-)2]"", where A -4.756, rg - 5.5 pixels, and a - 1.6 JKirsch at allbooeh 
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derived from the best-fit modeo Since Cash-statistics were developed for the modeling of 
unbinned data, we also modeled the unbinned spectrum for comparison. 

Two simple models were applied to the spectra. Model A is an absorbed power-law 
distribution. Model B is the same as A, except for the inclusion of a Gaussian component 
to model the Fe K emission at 6-7 keV (rest frame). Table [1] lists the best-fit parameters of 
each model and Figure [3] shows each set of spectra with the best-fit models. The significant 
improvement in fitting statistics of model B over model A suggests that there is indeed an 
emission line at an energy consistent with Fe Ka emission. However, since the number of 
counts is relatively low (especially when the data is binned to only 3 or 5 counts bin~^), 
the F-test cannot be used to determine whether the addition of the Gaussian component to 
model A is significant. The likelihood of the line being a result of statistical variations was 
tested using simulations. To this end, 10000 spectra were created using the fakeit command 
in XSPEC for each set of binned or unbinned data. The simulated spectra were created 
using model A. Then these spectra were fitted by both models A and B. If the line is a result 
of statistical variations, then one would expect a large fraction of the simulated spectra to 
be well described by model B. The fitting statistics were used to calculate Ac-stat(A-B) [or 
Ax^(A-B) for the 15 counts bin~^ data] which was then compared with the values presented 
in Tabled For the 15 counts bin"^ data, 1000 of 10000 (10.0%) had Ax^ greater than 3.76. 
This implies that model B (the inclusion of the emission line) is significant at the 90.0% level 
(a 1.6-cr detection). Similarly, the simulations show that the line is significant at the 96.87% 
level (313 out of 10000; 2.2a) for the 5 counts bin~i data, at the 93.5% level (507 out of 
10000; 1.8a) for the 3 counts bin'^ data, and at the 94.0% level (608 out of 10000; 1.9a) for 
the unbinned data. From these simulations, the line is significant to at least the 90.0% level. 

The 3 counts bin~^ data also suggest that the iron line can be decomposed into two 
narrower emission lines with centroid energies at 6.3 (EW ~0.6 keV) and 6.7 keV (EW 
~0.4 keV) in the rest frame. These energies are consistent with emission arising from neutral 



"^Gaussian statistics apply to data binned to at least 15 counts bin^^ while Poisson statistics apply to 
the data binned to at least 3 or 5 counts bin^^ and unbinned data. Since the difference of two Gaussian 
distributions remains a Gaussian distribution, a background-subtracted spectrum binned to at least 15 counts 
bin~^ retains the properties of a Gaussian distribution and can be modeled normally. However, the same is 
not true for a Poisson distribution. Therefore, the background cannot be simply subtracted for data binned 
to at least 3 or 5 counts bin~^ and unbinned data and then modeled. One way of treating the background is 
to model the background spectrum separately and then add the background model to the continuum model 
when fitting the source spectrum. For this paper, the background is modeled using a simple, relatively flat 
power law (P ^^1.0). This treatment of the background is applied to all modeling of data binned to at least 
3 and 5 counts bin^^ and the unbinned spectrum. A representation of the background spectrum and model 
is shown in the bottom panel of Figure 3] 
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iron and Fe XXV, respectively. The fitting statistics of tlie double-line model to the unbinned 
data is only slightly better than that of the single-line model. The detection of these narrow 
lines in the Fe K complex is significant at only the ~60% level based on 10000 simulations 
of the unbinned data. Therefore, the detection of the doublet needs to be confirmed with 
data of higher spectral resolution and signal-to-noise ratio. 

Our modeling and simulations show that Cash-statistics give consistent results for the 
unbinned spectrum and the spectra binned to at least 3 and 5 counts bin~^. Since Cash 
statistics were designed for unbinned spectra, we will use only the unbinned spectrum in 
subsequent modeling. The iron line is most prominent in the data binned to at least 5 
counts bin~^, we will use the spectrum binned to at least 5 counts bin"^ as a visual and 
qualitative check for the model of the unbinned data. 



3.3.2. AGN + Starhurst Continuum Models 

Aside from models A and B mentioned above, we modeled the unbinned spectrum with 
slightly more complex models to account for the possibility that a starburst may coexist 
with the AGN in F04103— 2838. Cautioning against over-interpreting data with only modest 
signal-to-noise ratios, even these more "complex" models were kept as simple as possible. 

The first model (model C) is a combination of absorbed power-law and MEKAL spectra 
(with metallicity fixed at solar) representing the emission from the AGN and starburst, 
respectively. The second model (model D) is a combination of two absorbed power laws, 
with one power law representing the AGN and the other representing the high mass X-ray 
binaries (HMXBs) associated with the possible starburst in this object. Finally, a third 
model (model E) was a combination of the two above mentioned models: a power law for 
the AGN, a power law for the HMXBs, and a MEKAL model for the hot gas. For all of 
these models, a Gaussian with centroid energy between 6 and 7 keV was included to model 
the iron line. 

While all of these models give better fitting statistics than the simpler power law models, 
only model C is a realistic fit to the data. Models D and E are rejected on the grounds that 
the best-fit power law values are physically unrealistic descriptions of AGNs. Therefore, we 
adopt model C as the "best-fit model" (Figure S]) and list the fitting parameters in Table [U 
This is perhaps not surprising given that ULIRGs are known from observations at optical and 
infrared wavelengths to show the presence of both an AGN and a starburst (e.g., Genzel et 
al. 1998; Kim, Veilleux, & Sanders 1998); F04103— 2838 does not appear to be an exception. 



Discussion 



4.1. The Soft Component 

The resuhs from the spectral fitting suggest that the soft X-ray (0.2-2 keV) flux is 
best described as thermal emission from hot gas with kT ~ 0.1 keV (T~1.2xl0^ K). This 
i s somewhat lower thari the range of gas temperatures (0.6-0.8 keV) found in LINERs 
(JGonzalez-Martin et al.l l2006r). The results for F04103-2838 is also somewhat lower than 



the results from iGrimes et al.l (120051 ) who performed a Chandra archival study of the soft 
X-ray emission from starburst galaxies ranging in luminosity from dwarf galaxies to ULIRGs. 
The authors found that the soft X-ray thermal emission of these starburst galaxies tends to 
fall in the temperature range kT ~0.25-0.8 keV with ULIRGs occupying the upper end of this 
temperature range. These large temperatures can all be attributed to powerful starbursts. 

The soft X-ray emission in F04103— 2838 is likely the result of thermal bremsstrahlung 
from a hot gas produced by the merger-induced starburst or by intrinsically extended soft 
X-ray emission heated by the AGN. If the ion density equals that of the electrons, the 
relationship between the electron density (ug) and luminosity of an emitting region of a 
given volume iV) is 

Lff ^ 1.7 X IQ-^^nlfV ergs s-\ (1) 

where / is the flUing factor for the hot gaqj- The non-AGN contribution of the nominal 0.2- 
2 keV luminosity from the best-flt model (model C) for F04103— 2838 is 1.6 x 10^^ ergs s^^. 
Assuming the emitting region is spherical with a diameter of <5'1j, the average electron 
density has a lower limit of ~0.19/~^'^ cm~^. This value is consistent with simulation results 
for the warm (10^'^ ^ T < 10^'^ K) component in the wind models of [Strickland fc Stevens 



fl2000[ ). 



Observationally, this hot gas compori e nt is difficult to probe because of its low den- 
sity and emissivity. [Strickland fc StevensI (J2000l ) performed hydrodynamic simulations of 
starburst- driven galactic winds with various ISM models. The authors found that, in gen- 
eral, the soft X-ray emission comes from gas with low fllling factors (10"'^ < / < 10~^; 
see also Cecil, Bland-Hawthorn, & Veilleux 2002; and Strickland et al. 2004a, 2004b for 



'^Equation [T] is based on equation (5.14b) and Figure 5.2 of lRvbicki fc LightmanI (jl979l ) for T 
the energy range of 0.2-2 keV. 



10^ K in 



^While the selection of a <5" emitting region is based on the spatial resolution of the telescope, it should 
be noted that the linear diameter of 5" at the distance of F 04103-2838 is less tha n a factor of two larger than 



the soft X-ray (0.5-2.5 keV) emitting region of NGC 6240 ( Komossa et al.ll2003l) . Therefore, the assumption 
of a <5" diameter is reasonable, even though it was chosen based on the instrument PSF. 
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observational constraints). Using these values for /, the electron d ensity of the h ot ga s in 
F04103— 2838 is ~0.6-5.9 cm"'^, consistent with values derived by iNetzer et al.l (120051 ) in 
NGC 6240. 

The soft X-ray emission detected in F04103-2838 may be thus the result of superwinds 
from the starbur st. X-ray superbub bles have been observed in Arp 220 (jlwasawa et al.ll2005l ) 
and NGC 6240 (INetzer et al.l 120051 ). Furthermore, powerful outflow events are now thought 
to take place in most ULIRGs (e.g., Rupke et al. 2002, 2005a, b, c, though their sample did 
not include F04103-2838). 



4.2. The Iron Feature 



F04103-2838 io ins the growing list of ULIRGs with Fe K detections 



2524 Jptak et al.l l2003h. Mrk 273 Jptak et al.ll2003[ ). and UGC 05101 Jimanishi et all 12003 



(llwasawa et al.ll2C)05r). Z11598-0112 (Paper I) F19254-7245 (IFranceschini et a 
2OO3I ), Mrk 231 (JMalonev fc Revnoldl IgOodJptak^et al.ll2003l : ISraito et alT 



e.g. 



Arp 220 



20031 : iBraito et al. 
F05189- 



200i). 



Ptak et al.ll2003l )]. supporting the view that an obscured AGN exists in many of these ob- 



jects. The presence of an AGN in F04103— 2838 was first suggested by Paper I based on the 
large hard X-ray to far-infrared flux ratio; the XMM detection of Fe K now indicates that 
the luminosity of this AGN has probably been underestimated. 



Few LINERs have detected Fe Ka lines. iTerashima et al.l (120021 ) studied a sample of 53 
LINERs and low-luminosity Seyfert galaxies using ASCA. Of the 21 LINERs in their sample, 
Fe emission lines were detected in only five galaxies (NGC 1052, NGC 3998, NGC 4261, 
NGC 4579, and NGC 4736). Of these five objects, only four (i.e. those excluding NGC 4261) 
have centroid line energies consistent within the uncertainties of the measurements with 
Fe Ka emission due to neutral iron (E ~ 6.4 keV). 



Three other LINERs have known Fe K detections; all t hree are power 



traluminous infra. r ed galaxies. These ga laxies are Arp 220 ( Iwasawa et al. 



ul lum inous or ul- 
2005I ). NGC 6240 



(iPtak et al.ll2003l : iKomossa et al.ll2003f ). and UGC 5101 (llmanishi et all l2003l : IPtak et al. 
2OO3I ). Chandra observations of Arp 220, the archetypal ULIRG, show an iron line at 



6.7±0.1 keV. Th is is consistent with e r nission due to Fe XX u p to Fe XXVI, but not neutral 
iron at 6.4 keV (llwasawa et al.l 120051 ). iKomossa et al.l (120031 ) detected Fe K emission from 
each of the two nuclei in NGC 6240. Their analysis showed that the iron lines in each nucleus 
are consistent with Fe Ka and Fe K/3 emissions. 

In Figure [5] we show the distribution of published Fe K equivalent widths of all LIN- 
ERs and ULIRGs known to have line emission. Arp 220, NGC 6240, and F04103-2838 
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appear to have iron emission with the greatest EW measurements of all the LINERs and 
ULIR Gs. These larg e Fe K features could be results of the blending of multiple narrower 
lines. iKomossa et al.l (120031 ) did not publis h the EWs o f the l ines from each of the nuclei in 
NGC 6240. The result quoted here is from iPtak et al.l (120031 ). The authors did not distin- 
guish Fe Ka emission from Fe K/3 emission and the EW measurement is likely dominated by 
the brighter southern nucleus alone. The large equivalent widths of the ULIRGs are telltale 
signs of obscured AGNs where line-of-sight columns of material exceeding lO^'* cm~^ prevent 
a direct view of the AGN; the 2-10 keV flux is dominated by light scattered off dust or 
electrons (e.g., Ghisellini et al. 1994; Krolik et al. 1994). The large amount of molecular gas 
(~ 10^ Mq pc-2) within 400 pc from the nuclei of NGC 6240 (e.g., Bryant & Scoville 1999) 
is sufficient to cause this obscuration. A similar explanation likely applies to F04103— 2838, 
although we are not aware of any CO measurements in this system. 

I nterestingly, the Fe K cor nplex in NG C 624(3 breaks up into a number of narrow lines. 
Both iNetzer et al.l (120051 ) and iBoUer et al.l (120031 ) detected Fe K lines due to neutral iron 
(^6.41±0.2 keV) , Fe X XV (6.68±0.02 keV), and Fe XXVI (7.01±0.04 keV) in NGC 6240. 
Komossa et al.l (l2003l ) also detected lines at 6.4 and 6.95 keV. The centroid energies of the 
lines due to neutral iron and Fe XXV in NGC 6240 are consistent with the respective cen- 
troid energies suggested by the doublet in the F04103— 2838 3 counts bin~^ data. Although 
simulations suggest the two line model is only significant at the ~60% level, a FWHM of 
~30000 km s^^ (a ~0.3 keV) seems too broad and the two component interpretation may 
be more likely. The Fe XXVI line in NGC 6240 is much fainter than the other lines so it is 
not surprising that we were unable to detect this feature in the modest signal-to-noise ratio 
data of F04103-2838. 

Despite their X-ray similarities, F04103— 2838 is ~2.5 times more infrared luminous than 
NGC 6240. These objects also differ in terms of IRAS /25//60 ratios (0.15 for NGC 6240 
and 0.30 for F04103— 2838) and merger state (NGC 6240 is in a pre-merger phase with a 
nuclear separation of ~ 1.3 kpc while F04103— 2838 is in the post-merger stage with a single 
coalesced nucleus). There is growing observational evidence (e.g., Veilleux et al. 2002, 2006; 
Ishida 2004; Dasyra et al. 2006a, b, 2007) and theoretical motivation (e.g., Hopkins et 
al. 2005) that mergers of gas-rich galaxies often produce "cool" (/25//60 < 0.2) luminous 
infrared galaxies that evolve into "warm" (/25//60 ^ 0.2) ULIRGs before becoming optical 
quasars. If this evolutionary sequence applies to NGC 6240 and F04103— 2838, the first 
object may actually be the precursor to the latter. 
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4.3. Energy Source of the ULIRG 

The lack of of short timescale variability (see §3.2p is to be expected if most of the 
primary X-ray flux is being absorbed or reprocessed. As discussed in §4.2, the large equivalent 
width of the iron line in F04103— 2838 implies the presence of a highly obscured AGN. 
It is very diflicult in such cases to estimate the intrinsic luminosity of the AGN without 
measurements of the > 10 keV flux from the buried AGN (e.g ., Mrk 231; Braito et al. 



2004). Here we follow the method of iMaloney fc Reynolds! (J2000l ) to estimate the intrinsic 
luminosity of F04103-2838. 

In their analysis of an ASCA observation of Mrk 231, they discussed two ways of es- 
timating the intrinsic AGN flux. The observed X-ray flux is due to a combination of two 



effects: reflection and scattering. iMaloney &: Reynolds! (J2000|) estimated the intrinsic AGN 
flux from the reflection and the scattering components separately. In their geometry, the 
observer has an obstructed view of the nucleus so the observed flux must be either scattered 
or reflected into the line of sight along which there is some amount of absorbing material. 
The reflected component is light from the central engine reflected off of the circumnuclear 
torus; the amount of reflection depends on the size of the reflecting surface and the compo- 
sition of the torus. On the other hand, the scattered component is light from the central 
engine (unobstructed by the torus) scattered into the line of sight. Based on their spectral 
fitting of the ASCA data, iMalonev fc Revnoldsl J2000h found that the X-ray flux of Mrk 231 
is scattering-dominated with 75% scattered and 25% reflected light. 

Due to the low si gnal-to-noise ratio of our d ata on F04103-2838, the same spectral 
fitting as performed by !Maloney &: Reynolds! (|2000! ) could not be done. The large equivalent 
width of the Fe Ka line (~1.6 keV) above 1 keV suggests a reflection-dominated spectrum. 
However, the width of the line implies it could be a blend of narrower Fe Ka and ionized 
iron emission lines (as suggested by the 3 counts bin~^ data). If this were the case, the 
Fe Ka EW may be more consistent with a scattering dominated spectrum. Therefore, we 
will consider two cases: (1) the majority of the observed flux is due to reflection and (2) the 
majority of the observed flux is due to scattering to estimate the intrinsic X-ray luminosity 
of the AGN. 

After correction for absorption, the nominal 0.2-10 keV flux of the buried AGN in 
F04103-2838 derived from our best-flt model (model C) is 1.83 x 10^^ ergs s~^. In the flrst 
scenario, we will assume the reflection component is 75% and the scattering component 
is 25% of the total observed flux. This implies that -Lsnattered = 0.45 x 10"^^ ergs s~^ an d 
i^refiected = 1-38 X 10^^ ergs s"^ for the AGN in F04103-2838. In lMalonev fc Revnoldsl J200oh . 



the luminosity from the reflected portion is scaled up by a factor of 25 in their modeling of 
the reflection process. The reflection process differs for different galaxies; it depends on the 
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ionization state of the mirro r and the steepness of the photon index of the central black hole. 
Maloney &: Reynolds! (120001 ) assumed reflection from neutral material, a reflecting fraction 
of 10%, and th e canonical value of the pho ton index due to an AGN (T = 1.8). The scaling 
factor used by iMaloney &: Reynolds! (120001 ) corrects for the flattening of a spectrum with T 
= 1.8 to r ~1.1 (based on a single absorbed power-law model) for Mrk 231 due to reflection. 
The correction factor of 25, therefore, is a maximum correction factor. The minimum scaling 
factor is 10 to simply correct for a reflecting surface fraction of 10%. We will conservatively 
assume this minimum scaling factor of 10 for the reflectio n component. For the scatt ering 
component, we will assume the same scattering fraction as !Maloney &: Reynolds! (|2000! ) (i.e. 
1% for electron scattering). After the corrections, the intrinsic 0.2-10 keV luminosity of the 
AGN in this scenario, where reflection dominates the observed flux, is 1.4 x 10^^ ergs s~^ 
(from reflection) and 4.5 x 10^^ ergs s~^ (from scattering). Therefore, the total reflection- 
and scattering-corrected luminosity in the 0.2-10 keV band is 5.9 x 10^^ ergs s~^ if we assume 
reflected light dominates the observed spectrum. 

Similarly, for the second scenario where the majority of the observed flux was scattered 
into the line of sight, we will assume the reflection component is 25% and the scattering 
component is 75% of the total observed flux. This implies that the intrinsic 0.2-10 keV 
luminosity of the AGN is 4.5 x 10^^ ergs s~^ (from reflection) and 1.4 x lO'*'^ ergs s~^ (from 
scattering). Hence, the corrected 0.2-10 keV luminosity is 1.4 x 10^"^ ergs s~^. It should be 
noted that in both cases we considered, the luminosity from the scattering portion dominated 
the total after the corrections. 



Thus, with the assumptions made above, the intrinsic 0.2-10 keV luminosity of the 
AGN ranges from 0.6 to 1.4 x 10^^ ergs s~^. This range in luminosity overlaps with that of 
quasars (~10^'^ ergs s~^; e.g., Elvis et al. 1994; Piconcelli et al. 2005) and is similar to that 
of NGC 6240 (~ 0.7 — 2 x 10^^ ergs s~^, after correction for an HI column density of 1 — 2 x 
10^4 cm-2; Vignati et al. 1999). The ratio log[L2_io kcv/Lm] for F04103-2838 corrected for 
scattering and reflection is -2.2 to -1.7. These values fall precisely within the range found in 
radio-quiet PG quasars (-3 to -1; Sanders et al. 1989). Ass uming F04103-28 38 has the same 
X-ray to bolometric luminosity ratio as radio-quiet QSOs (JElvis et al.l!l994l . Lx/Lboi ~3%), 
the AGN contribution to the bolometric luminosity of F04103-2838 is ~15-38%. Therefore, 
within the large uncertainties, the AGN in F04103-2838 does not dominate the total energy 
output of the galaxy. 
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5. Summary 

The results from our analysis of the XMM-Newton spectrum of the 1-Jy ULIRG/LINER 
F04103— 2838 can be summarized as follows: 

1. The soft (0.2-2 keV) X-ray flux of F04103-2838 is attributed to hot gas with kT 
~0.1 keV. This temperature is similar to that derived in other starburst galaxies and 
LINERs. The electron density in F04103-2838 is ~0.6-5.9 cm^^, consistent with the- 
oretical predictions and observational estimates in wind systems. 

2. An Fe Ka line located at ~ 6.4 keV with an equivalent width of ~1.6 keV is detected 
in F04103— 2838. The line could be intrinsically broad or could be made up of two 
narrow lines located at rest frame energies of ~ 6.3 and 6.7 keV but this decomposition 
is only significant at the ~60% level, so it needs to be verified with higher resolution 
spectra. 

3. The large equivalent width of the Fe Ka line suggests that the AGN is Compton-thick. 
Using simple assumptions, we estimate that the intrinsic 0.2-10 keV luminosity of this 
AGN is 0.6 - 1.4 xlO^^ ergs s^^. If these assumptions are correct and the galaxy has a 
QSO-like X-ray to bolometric luminosity ratio, the AGN detected by our observations 
does not dominate the bolometric luminosity of F04103— 2838. 

4. The X-ray spectral characteristics of F04103— 2838 are strikingly similar to those of the 
local luminous infrared galaxy NGC 6240. Given the similarities in X-ray properties but 
differences in merger state and in infrared color and luminosity, objects like NGC 6240 
could conceivably be the precursors of ULIRGs like F04103— 2838. 
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^Model A: Absorptioncaiactic x Absorptiongourco x PL. Model B: Absorptioncaiactic x Absorptiongource x (PL + Line). Model C: 
AbsorptionGaiactic X [MEKAL + Absorptionsourco x (PL + Line)], -where MEKAL is the Me-we, Kaastra, & Liedahl thermal plasma 
model (see the XSPEC manual for details), PL is a po-wer-law model representing the AGN, Line is the Fe K emission line -with 
a Gaussian profile, AbsorptionGaiactic is the absorption from Nh, Galactic — 2.45 x 10^° atoms cm^^, and Absorptionsourco is the 
intrinsic absorption within the source. 

^Intrinsic (i.e. "within the galaxy) column density in units of 10^^ atoms cm~^. 

'^Fe K line energy (rest frame), "width, equivalent "width, and thermal gas temperature all in keV. 

■^Fitting statistics per degrees of freedom. Gash statistics are used for unbinned spectra and spectra binned to at least 3 and 5 
counts bin"'^, "while x^ statistics are used for spectra binned to at least 15 counts bin~^. 

'^ Absorption-corrected flux in units of 10^^'* ergs cm^^ s^^. The AGN value includes the flux from both the po"wer-law component 
and the iron line. 



''Absorption-corrected luminosity in units of 10^^ ergs s ^ . The AGN value includes the flux from both the power-law component 
and the iron line. 
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Fig. 1.— Lejt: EPIC mosaic image of F04103-2838 smoothed with a Gaussian fF WHM~5^0 



and d isplayed on a hnear grey scale. The contours are optical R-band data from iKim et al. 



fl2002[ ). Right Comparison of unsmoothed EPIC PN PSF at 1 keV with observed 0.2-2 keV 
radial profile of F04103-2838. The dashed line represents the theoretical PSF while the solid 
line represents the PSF broadened due to uncertainties in the correction for pointing drift of 
the telescope. The absolute pointing drift (APD) error is conservatively assumed to be 3", 
the u pper lirn i t (see XMM-Newton Observer's Handbook). The error bars were calculated 
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Gehreld ( 1l986l ). The source is unresolved within the uncertainty of the measurements. 
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Fig. 2.— The 0.2-10 keV and 2-10 keV light curves for F04103-2838. The sohd crosses 
denote the source count rate, while the dotted crosses denote the background count rate. 
Background-subtraction was not applied to the source spectrum. The time bins are each 
5234 seconds. The error bars are Poissonian counting errors calculated following iGehrels 
(1l986l ) at the 84% significance level. Within the errors, the source is not variable on the 
~20 ksec time scale of our observation. The lack of short timescale variations is expected of 
sources where most of the primary X-ray flux is absorbed or reprocessed. 
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Fig. 3. — EPIC PN spectrum and best-fit models to F04103— 2838 with different binnings: 
> 15 counts bin~^ (top panels), > 5 counts bin~^ (middle), and > 3 counts bin~^ (bottom). 
The unbinned spectrum was modeled but is not shown here. Model A (left panels) is a 
simple absorbed power-law distribution; model B (right) is the same as A, but includes a 
Gaussian component to model the Fe K emission. The best-fit model parameters are listed 
in Table [H The iron line is most prominent in the data binned to 5 counts bin~^ and the 
Fe K doublet may be present in the 3 counts bin~^ data. X-axis of the figures represents 
energy in the observer's frame. 
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Fig. 4. — EPIC PN source and background spectra of F04103— 2838 binned to at least 5 
counts bin~^ with the best-fit unbinned model (C) applied. The X-axis of the figures repre- 
sents energy in the observer's frame. Top: The source model includes a MEKAL component 
for the thermal emission, a power- law component to represent the AGN component, a Gaus- 
sian component to model the Fe K emission, and a relatively flat power-law (F ~1.0) for the 
background. The model parameters are listed in Table [TJ Bottom: A binned background 
spectrum with the background model used in the modeling of the unbinned spectrum and 
spectra binned to at least 5 and 3 counts bin~^. No significant features are seen in the 
background spectrum. 
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Fig. 5. — Distribution of equivalent widths of Fe K emission features in LINERs and 
(U)LIRGs versus absorption-corrected 2-10 keV to far-infrared flux ratio. In general, the 



(U)LIRGs have higher EWs than IR-fain t LINERs. O 



was determined to be Compton-thick 



figure are drawn f rom iTerashima et al. 



by. 



Ptak et all (l2003f ). lBraito et al.l fl2004 ). 



Braito et al. 



J2002h. 



these, F19254-7245 (EW~2 keV) 
fl2003r) . The values included in this 



Braito et al. 



fl2003h. 



Imanishi fc Terashimal (120041 ) 



manishi et al. 



Iwasawa et al.l (120051 ) 



( 2003 ) 



Paper I, and this work. The equivalent widths for the four LINER (U)LIRGs included in 
this sample ar e labeled. The F04103 -2838 value plotted is derived from our best-fit model 
(C). Note that JKomossa et al.l (120031 ) detected Fe K emission f rom e ach of the two nuclei in 
NGC 6240. Therefore, the single value quoted by iPtak et al.l (120031 ) is a sum of the Fe Ka 
and Fe Kj3 emission due to neutral iron, likely dominated by the brighter southern nucleus. 
The value quoted for Arp 220 may also be due to a blend of emission lines arising from 
ionized iron (Fe XX up to Fe XXVI). The value (~0.5 keV) quoted for the ULIRG Mrk 463 
is a sum of the emission due to neutral iron and Fe XXV. 



